Large N Scaling Behavior of the Lipkin-Meshkov-Glick Model 
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We introduce a novel semiclassical approach to the Lipkin model. In this way the well-known 
phase transition arising at the critical value of the coupling is intuitively understood. New results 
- showing for strong couplings the existence of a threshold energy which separates deformed from 
undeformed states as well as the divergence of the density of states at the threshold energy - are 
explained straightforwardly and in quantitative terms by the appearance of a double well structure 
in a classical system corresponding to the Lipkin model. Previously unnoticed features of the 
eigenstates near the threshold energy are also predicted and found to hold. 
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The Lipkin-Meshkov-Glick (LMG) model, originally 
introduced in nuclear physics has found applications in 
a broad range of other topics: statistical mechanics of 
quantum spin systems [2( , Bose-Einstein condensates Q 
as well as quantum entanglement Q , to name but a few. 
The continued interest in this system arises from the fact 
that it is an exactly solvable 0.0 many-body interacting 
quantum system as well as one of the simplest to show a 
quantum transition in the regime of strong coupling. 

This transition is by now well-understood: the ground 
state becomes degenerate and a macroscopic change in 
the ground state energy takes place. Furthermore, at the 
transition value of the coupling, the density of states di- 
verges at the ground state energy as the number N of 
interacting particles becomes large. The energy gap was 
recently found to vanish as N^ 1 ^ 3 . Furthermore, in 
this model a novel type of phase transition has recently 
been discovered for strong values of the coupling param- 
eter. Indeed, in this regime the spectrum is divided by 
a critical energy E c , where the behavior characteristic 
of strong coupling holds below E c , while above E c the 
system reverts to the kind of behavior found below the 
phase transition. At the critical energy, the density of 
states is again found to diverge as the particle number N 
goes to infinity. This divergence has been conjectured to 
be of the logarithmic type. 

In recent papers 
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different approaches have been 
used to explain these results. The continuous unitary 
transformation technique (flow equations) was ap- 
plied to obtain the spectrum reliably for large N values. 
An investigation of the singularities of the spectrum lal 
in the complex A-plane (the exceptional points 
sheds more light upon the complexity of the limit prob- 
lem but a final answer about the limit attained has not 
been given. These singularities have been recognized as 
an essential mechanism to invoke the phase transition 
in a seminal paper by Lee and Yang [12J and their sig- 
nificance for the partition function is discussed more re- 
cently again in 0, ^| . A different approach starts 



with the bosonization method using higher orders of pow- 
ers in 1/N of the Holstein-Primakoff representation and 
then applies the flow equation technique. New results 
have been obtained in this way, in particular the correct 
analytic behavior of the level distance as a function of TV 
at the critical point. Earlier attempts [15j have revealed 
different results such as the form of the wave functions 
beyond the phase transition. 

In this Letter we introduce a novel semiclassical ap- 
proach to the LMG model. It readily explains the above 
features, determines the precise value of the critical en- 
ergy as a function of the coupling, proves the logarithmic 
divergence of the density of states near the critical energy 
as well as successfully predicts certain previously unno- 
ticed behavior of the eigenstates near the critical energy. 
Finally, we obtain a qualitative understanding of all es- 
sentials of the model: the classical model we introduce 
has a double well structure above the phase transition, 
and the critical energy can then be identified with the 
separatrix energy. The approach given here also shows 
easily both the nature of the phase transition as a func- 
tion of the coupling parameter as well as the scaling with 
N of the vanishing gap at the critical coupling, which 
was previously shown |j| to scale as TV -1 / 3 . 

We here recapitulate the basics of the model and dis- 
cuss the essential properties for large values of N. It is 
given in terms of 2j + 1 = N + 1-dimensional represen- 
tations of the SU(2) operators Jk, k follows 

H(\) = J X + ±(J*-J$). (1) 

Here the interaction is scaled by N to ensure that H is 
extensive. In this form the model has a phase transition 
just beyond A = 1, the larger N the closer the transition 
point at A = 1. This has been discussed under various 
points of view in the literature, see e.g. 0E|. 

The Hamiltonian allows reduction into two spaces: m 
integer and m half-integer, with m the eigenvalues of J 2 ; 



2 



it corresponds to N even and odd respectively and is de- 
noted as parity. For A > the even and odd levels are 
obviously separated and remain so for all A < 1 while 
the levels become degenerate (up to terms vanishing ex- 
ponentially fast in N) for A > 1. The phase at A < 1 
is called the normal phase, while the symmetry (parity) 
breaking phase at A > 1 is called the deformed phase. 
Recent calculations |.9j apply non-perturbative flow equa- 
tions allowing to obtain the spectrum for arbitrarily high 
yet finite values of N. These have established || the exis- 
tence of the phase transition in energy referred to above: 
the states having energy below a certain threshold be- 
have as states of the deformed phase, whereas higher in 
the spectrum the states become undeformed again. 

Since the commutator [H(\), (J) 2 ] vanishes, we confine 
ourselves to a fixed value of j = N/2. For large N we 
consider the Hamiltonian (JTJ on the sphere of radius j = 
N/2. In other words, we rewrite the Hamilton operator 
UJ as a classical Hamilton function 

H = y ( - sin 6 cos </> - ~ (cos 2 9 - sin 2 9 sin 2 <f>)j (2) 

having introduced the polar angles as 

N 

J z = — — sin 9 cos <p 
N 

J x = — cos 9 
2 

N 

J v = — sin 9 sin 6. 



Note that the transition to a classical Hamiltonian has 
also been achieved in a different way using the coherent 
state approach developed in [3| ■ At this point we notice 
that, with n — cos 9, the Poisson bracket 



(3) 



suggests how to quantize the Hamilton function of the 
single particle problem in the two canonical conjugate 
coordinates <j) and li. It can be written as 
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2H 



A 



A1 2 cos ( /»--(^-(l-^)sin^). (4) 



To obtain information about the ground state and the 
low lying states we expand 2H/N around its minimum 
which is found at 



sin <po = 



Mo = 



(A < 1) 

±Vl- A- 2 (A > 1) 



and the corresponding minimum values of H at 
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Expanding around the minimum reveals in fact all essen- 
tial features. Keeping to lowest order terms we obtain 
around fi = <f> = 



K = -1 



1 + A 



1 - A 



(8) 



In this form, quantization is straightforward. Based on 
P|) we identify <f> and li with momentum and position, re- 
spectively, and 2/N with h, i.e. we use the usual canonical 
commutation relations for li and <f>. With this identifica- 
tion the Hamiltonian (JHJ represents a quartic oscillator 
which behaves for A < 1 basically like a traditional os- 
cillator with a harmonic spectrum (for the lower states) 

Ek ~ fc\/l — A 2 , fc = 1,2, For A > 1, the lower states 

must be determined at the minimum in li around the val- 
ues iio as given in ©. With e = li — (j,q the expansion 
of the Hamiltonian K yields, up to the constant given in 
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For low- lying levels, the harmonic spectrum is again ob- 
tained from © with frequency y^(A 2 — 1) well-known 
from previous work [3, Il9| . 

To evaluate the average energy density, the features of 
which have been the object of recent work we exploit, 
as we are in the semi-classical regime, the WKB relation 



S(2E k /N) = 2tt [k 



(10) 



Here S(E) denotes the action corresponding to the 
Hamiltonian K . By differentiation we obtain 

■whN 2?r . , 

AE = Eu+x -E k = = = 11) 

+ T(2E/N) T(E) V ; 

where T(E) is the period of the orbit with respect to K 
as a function of energy and E is (Ek+i + Ek)/2. Use is 
made of T(E) = dS/dE. 

The last relation (|llfl contains virtually all basic infor- 
mation of the Lipkin model. Firstly, for A < 1 nothing 
dramatic happens: the density of states merely changes 
smoothly as the period varies. An important fact should 
still be noted: the bosonization approach predicts, a con- 
stant AE in this case. We see here that this does not hold 
over the entire energy range: the period varies smoothly 
over large energy scales, and so does the average energy 
spacing; it increases with the energy. 

Secondly, for A > 1 there is an energy, in fact the 
separatrix, being situated for K at the value —1, where 
the period T(E) diverges. To estimate the average energy 
spacing in this region, we consider trajectories near the 
classical separatrix where they spend a long time near 
the unstable equilibrium, that is, in a region in which (j) 
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FIG. 1: Asymptotic behavior of the distance of levels of one 
parity at the transition point as a function of A. The points 
are the numerical fits for 500 < N < 1500 when fitted to 
/(A)/ln(A); the solid curve is 2ttVA 2 - 1. 



as well as fi are small. One has for K approximately 



1 + A 



A- 1 



(12) 



Since this dominates the divergence of T(E) near E = 
— 1, one finds that T(E) is approximately given by 
\n\2E/N + 1|/VA 2 - 1. But 2E/N + 1 is itself of the 
order 2AE/N and hence of first order in 1/AT. We thus 
obtain 



AE = 



2ttVA 2 - 1 
inJV 



(13) 



This corresponds to the high density of states observed 
in Q for specific values of A > 1. It was found that 
a change occurred between the two regimes below and 
above a certain A-dependent energy E C (X): for low ener- 
gies E < E C (X), the states were deformed, the order pa- 
rameter was non-zero and an odd-even degeneracy was 
observed. For E > E C (X), all these phenomena disap- 
peared and a normal regime, similar to A < 1 was recov- 
ered. The transition region between those two regimes 
had the typical signature of high density of states. From 
(|ll|l it becomes clear that the lower portion of the two 
regimes correspond to bounded motion in one well break- 
ing parity symmetry. We mention that the tunneling 
between the left hand and right hand wells determines 
the splitting to be ~ exp(— const jK) ~ exp(— const N). 
For higher energies, the corresponding classical motion 
is above the wells and symmetry is restored. The two 
different regimes are separated by the separatrix with its 
high density of states. Formula (|13fl has been verified 
numerically as shown in Fig. 1. The region of high den- 
sity is predicted to occur at energy — 1 corresponding to 
energies above the ground state around (A + A -1 — 2)/2 
which again is verified numerically. 

As the third major result, i|ll|) provides the leading 
analytic behavior of the spectrum at the transition point 



(A = 1). Since the Hamiltonian is quartic for (classically) 
low energies, one can evaluate AE using the formula for 
T(E) for a quartic oscillator valid for energies high up 
in the spectrum but still small with respect to AT. It is 
given by T(E) = const. • iT 1 / 4 , fr om which 



AE = const. • (E/N) 



1/4 



follows, and it entails 
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The AT -1 / 3 behavior has been obtained recently and 
the fc 4 / 3 behavior is confirmed numerically (see Fig. 2). 
We stress the non-uniform nature of the limit played by 
the critical point A = 1 when comparing (|14|l and (|13|) . 

It is obvious that the spectra of the two Hamiltonians 
JU and J2J) differ. This is due to the issue of ordering: 
in order to make sense of (J2Jl, we must specify in which 
way we order \i and <p to obtain a self-adjoint operator. 
There is no unique prescription for this so that an un- 
known difference exists between the two Hamiltonians. 
It is, however, known that if care is taken, these errors 
are of order h 2 , that is, of order N~ 2 . From this follows 
that we also expect the singularities to be different. We 
recall that they are associated with the critical point and 
the transitions for A > 1. As the semi-classical treatment 
preserves these basic features it is expected that at least 
the qualitative pattern of the exceptional points remains. 
In fact, there is (i) the special feature at A = 1 (an ac- 
cumulation point for N — > co (ii) a high density of 
EPs near to the separatrix 0] , that is for energies around 
E C (X), (iii) the absence of singularities near to real values 
for A < 1, and for A > 1 for energies sufficiently distant 
from E C (X). 

Of interest in the semi-classical treatment is the be- 
havior of the wave function at the phase transition for 
A > 1. In accordance with the long dwelling time clas- 
sically at the saddle point, there is the phenomenon of 
super-scarring for the wave function. In fact, this has 
been shown generally [2(J for the occurrence of such dou- 
ble wells. It arises for the specific values of k where the 
separatrix itself satisfies the WKB condition Hl()(l . The 
wave function then shows a dramatic concentration at 
the saddle-point of the Hamiltonian (or at the maximum 
in configuration space): in particular, it can be shown 
that there is an interval, the length of which goes to zero 
with h, that is with 1/AT, in which the whole wave func- 
tion is asymptotically concentrated. This is reflected in 
the eigenvectors of Q associated with the eigenvalues 
at the minimal gap: in the basis of J z they become rel- 
atively concentrated in the first few components, that 
is around fi — 0; while the number of substantial com- 
ponents increases with N, the first twenty components 
exhaust the norm by about 50% irrespective of N. This 
has been tested for A = 1.1, 1.5 and 2.0 where the tran- 
sition occurs at k w AT/120, AT/16 and AT/8, respectively. 
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FIG. 2: Log- log-plot of Ek versus k. For clarity only every 
eighth point of the first 500 levels are taken (TV = 5000). The 
straigth line fits the slope 4/3 with 1% accuracy. 



scarring at the saddle point in phase space are presented. 
The qualitative behavior of the singularities of the spec- 
trum (the exceptional points) seems to be preserved in 
the semi-classical approach. 
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The effect is rather significant - the more so the nearer 
A to unity - in that the first twenty components of the 
neighboring wave functions contribute appreciably less 
to the total norm. In principle, wave functions can be 
directly determined semi-classically. As this is not the 
major focus of this paper we only outline the procedure. 
Switching to the usual polar coordinates defined around 
the z-axis but keeping the definition of fj, as before, the 
renormalized Hamiltonian K reads 

K = fi+ ^(1- /r 2 ) cos 20. (16) 

In these variables, the usual spherical harmonics are 
eigenf unctions of the operator p, = (2/iN)d/d<j> with an 
eigenvalue proportional to that of J z . The semi-classical 
eigenf unctions of (|16|l in the eigenbasis of fi are then ex- 
pressed as a function of the energy by means of stan- 
dard WKB formulae for the one-dimensional Hamilto- 
nian (|16|l . The expressions fail, of course, in the usual 
manner near the turning points of i|16fl . 

To summarize: using the semi-classical version of the 
original model Q leads to J2J.Th.c- WKB approximation 
then yields the expansion JSJ . Most information can then 
be extracted from Ijlll) being based on H10|) . There is (i) 
the qualitative result about the phase transition occur- 
ring for A > 1 as discussed in |q. This includes the 
deviation from the strict equidistant level sequence for 
large N and for A < 1 as well as the exponential (in N) 
separation of the degenerate levels for A > 1. There is 
(ii) the expression l|13(l for the level distance at the tran- 
sition point for A > 1; there is (iii) the finding (|15|l at 
A = 1. The apparent contrast of results (ii) and (iii) un- 
derlines once again the nonuniform nature of the large N 
limit at A = 1. Note that all these results have been 
confirmed numerically. Additional results referring to 
the semi-classical wave function, in particular the super- 
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